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By John DeYoung 


SUMMARY 

The boundary conditions for low-aspect-ratio wings are approxi- 
mately satisfied by a two-dimensional vortex sheet that extendB to 
infinity fore and aft. Since the downwash is constant fore and aft, 
the boundary condition that no flow passes through the wing or flat 
plate is satisfied simply by the familiar downwash integral used in 
line— theory loading methods. A theory is developed herein, based on 
the foregoing principle which (for small aerodynamic angles and for 
plan forms that have all points of the trailing edge at or behind the 
/ upstream line of maximum span) will predict the spanwise loading due to 

arbitrary spanwise angle-of-attack distribution, including (l) uniform 
angle of attack, (2) that induced by rolling, (3) Bpanwise variation of 
* twist or camber, and (4) that due to deflection of arbitrary control 

surfaces of low aspect ratio. The spanwise loadings are independent of 
the plan forms considered. 

It is found for plan forms having unswept trailing edges that the 
spanwise loading distribution depends only on the slope of the region 
at the wing trailing edge. Thus, in the limit case of zero aspect ratio 
of wings and control surfaces, the spanwise loading is independent of the 
position of the hinge line along the wing chord. In addition, the total 
loading due to any chordwise camber distribution depends only on the 
slope at the wing trailing edge, provided the curvature at the trailing 
edge is small as compared with the slope. 

From considerations of the linearized— compressibleHTlow equation, 
the aerodynamic values given by zero-aspect— ratio theory are also those 
of moderate-aspect-ratio wings at the speed of sound. The results are 
applicable for low values of the product of the Prandtl-Glauert com- 
pressibility parameter and aspect ratio. 


INTRODUCTION 

The aerodynamic characteristics of pointed-wing plan forms with 
aspect ratio approaching zero have been considered by R. T. Jones 
(reference l) and extended by H. S. Ribner (reference 2) to include all 
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the stability derivatives for the pointed triangular wing. The triangu- 
lar wing of references 1 and 2 is treated on .the assumption that the 
flow potential in planes at right angles to the long axis of the airfoil 
is similar to the corresponding two-dimensional potential. This theory 
gives the surface loading for the flat plate and can he extended to the 
computation of all the important aerodynamic characteristics. 

Reference 3 obtains surface loadings for swept wings with the 
linearized-compressible— flow equation of supersonic theory, which, at 
the speed of sound, results in equivalent plan forms that have zero 
aspect ratio. Reference 3 gives the lift coefficient and induced drag 
for a range of constant-chord swept wings, which results are applicable 
for low-aspect-ratio swept wings throughout the speed range and moderate— 
aspect-ratio wings at the speed of sound. 

A theory for the rectangular wing of low aspect ratio has been given 
by W. Bollay (reference 4 ). This work considers the flow leaving the 
side edges of the plate at all angles of attack and thus approximates the 
early separation occurring on blunt plan forms. The results indicate 
that lift varies nonlinearly with angle of attack as aspect ratio becomes 
small. It can be assumed by the theory of reference 1 that the flow 
sticks to the plate, an assumption which is more valid for pointed than 
for blunt plan forms since loading is continually generated on the swept 
leading edge. Thus, lift varies linearly with angle of attack for a 
greater range of angle of attack for the pointed than for the blunt plan 
form. 


There has been, up to the present, no published work presenting 
control— surface and twist effects for low-espect-ratio wings or for 
moderate-aspect— ratio wings at sonic speeds. The primary purpose of the 
present report is to extend the theory of reference 1 to provide a means 
of predicting aerody n am i c characteristics resulting from the deflection 
of arbitrary control surfaces of low aspect ratio on low-aspect-ratio 
wings. Since the linearized-compressible— flow equation indicates that 
low values of the product of the Prandt 1-Glauert parameter and aspect 
ratio are satisfied by low-aspect-ratio theory, the results will be 
applicable throughout the speed range for low-aspect— ratio wings and at 
sonic speeds for moderate-aspect-ratio wings. 


SYMBOIS 


A aspect ratio (b 2 /s) 

a n coefficients of the trigonometric series for spanwise loading 

b wing span measured perpendicular’ to the plane of symmetry, feet 
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wing chord, feet 1 


per radian 


lift due to control— surface deflection 


ro lling moment due to rolling 


per radian 


lift due to angle of attack f > P® r radian. 

/dCr \ 

lift due to control— surface deflection ( 1, per radial 

r ac, i 

ro lling moment due to rolling I ^pb/gy) ’ J * per rai ^ an 

/a Ci\ 

ro lling moment due to control-surface deflection l ■ ) , 


per radian 


per radian 


wing lift coefficient 


•m il ing -moment coefficient 


local lift coefficient 


( total lift N 

\ qS J 

lent ( ro ^- 1rig ^^nt \ 

611 \ qSb / 

/ local lift ^ 


induced drag coefficient 

spanwise loading coefficient 

/ T \ ' ' y 

circulation i — ) 

\bV / 


or dimensionless 


circulation 


Mach number 


integer coefficients of loading series 
rate of rolling, radians per second 


rolling helix angle, radians 


q. free-stream dynamic pressure, pounds per square foot 

S wing area, square feet 

V free— stream velocity, feet per second 


•‘•Quantities measured parallel to the plane of symmetry. 


k 
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v induced velocity, normal to the wing, positive for downwash, 

feet per second 

y lateral coordinate pertaining to loading, measured from wing 

root perpendicular to plane of symmetry, feet 

y lateral coordinate pertaining to downwash, feet 

a constant spanwise angle of attack, radians 1 

a(cp) variable spanwise angle of attack, radians 1 

wing root angle of attack when total lift equals zero, 

radians i 

0 compressibility parameter 

T spanwise circulation, feet squared per second 

81,82 arbitrary control-surface angles for each semispan, radians 1 
e(cp) wing twist relative to root chord, radians 1 
i) dimensionless lateral coordinate pertaining to loading ( -^7- 


dimensionless lateral coordinate pertaining to downwash 


(*) 


^c.p. semispanwise center^-of-pressure location 


parameter defining position of control-surface spanwise edge, 
radians 


taper ratio 


tip chord 


\ root chord / 

sweep-angle parameter, positive for sweepback ^tan - 1 , 

degrees 


Subscripts 

integer pertaining to span station 


1 See footnote 1, page 3 
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a ailerons 

f flaps 


THEORY FOR W33R5S OF ZERO ASPECT RATIO 
Plait-Form Range and Boundary Conditions 

The surface loading on a wing can he represented hy a vortex sheet. 
Hie vortex— sheet distribution over the surface of the low-aspect-ratio 
triangular wing is pictured in figure 1(a). As the aspect ratio 
approaches zero the distribution of figure 1(a) becomes the simple two- 
dimensional fore and aft vortex distribution pictured in figure 1(b). 
Various wings having unswept trailing edges and approaching zero aspect 
ratio are pictured in figure 2. The vortex distributions of these wings 
are (s imilar to the wing of figure 1(a)) increasingly approximated by 
the two-dimensional vortex distribution of figure 1(b) as aspect ratio 
approaches zero. Hence, as aspect ratio approaches zero, the dowrrwash 
f normal to the wing, due to the vortex sheet of the wing, approaches 

that due to a two-dimensional vortex sheet. The plan-form boundary ■ 
conditions that there be no flow through the plate can be easily satis— 

- fied by the two-dimensional vortex sheet because the downwash in the 

longitudinal direction is constant. The spanwise loading distribution 
is equal to the vortex strength distribution in the sheet. 

From the linearized— compressible-flow equation, the aerodynamic 
characteristics of wings approach those given by zero aspect ratio as 
the wing nears the speed of sound. The compressibility equation 
requires that the longitudinal coordinate of a point in the plan form 
effectively increase as l/f3. With this condition (referred to as the 
Prandtl transformation), wings approach zero aspect ratio in the same 
manner as those depicted in figure 2. Thus, at the speed of sou nd , plan 
forms having unswept trail ing edges and moderate aspect ratios will also 
have their boundary conditions satisfied by the two-dimensional vortex 
sheet . 

The two-dimensional vortex sheet can satisfy boundary conditions on 
the wing downstream from any chordwise angle— of-attack changes such as 
those due to control surfaces and camber. In the case of a localized 
change in angle of attack, the loading will actually be developed in a 
region near the bend in the wing surface, and in this region the 
boundary condition will not be satisfied by the simple two-dimensional 
vortex sheet. However, all loading changes that occur in any region of 
the wing are taken into account by the strength of the trailing vortex 
sheet, which must take up the configuration to satisfy the boundary 
» conditions along the relatively infinite wing chord downstream. Thus, 

for wings of zero aspect ratio having unswept trailing edge, the loading 
depends only on the final angle of 'attack attained by the wing in the 
region of the trailing edge, provided the curvature of the wing chord is 
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s mall in ratio to the slope. The spanwise loading distribution is then 
independent of the chordwise location or distribution of a control- 
surface hinge line in the limiting case of zero aspect ratio or for 
moderate aspect ratio at the speed of sound. Also, for such wings, the 
total wing loading due to any chordwise distribution of camber depends 
only on the final angle of attack of the wing in the region of the 
trailing edge, provided the ratio of curvature to slope is small in this 
region. 

Plan forms that approach zero aspect ratio by the Prandtl trans- 
formation, having other than an unswept trailing edge, can have the 
boundary condition satisfied by the two-dimensional vortex sheet, pro- 
vided all points of the trailing edge remain at or behind the upstream 
line of ww-TiTmim wing span. However, for these plan forms, the two- 
dimensional vortex sheet will only satisfy boundary conditions of 
chordwise changes in angle of attack that are at or ahead of a span- 
wise line through the max-timim upstream point of the trailing edge. The 
reason for the limitation of wing plan forms is that the two-dimensional 
vortex sheet is applied to the full span of the wing and, if a t railing- 
edge indentation exists ahead of the upstream line of maximum span, part f 

of the wing loading would be off the wing. The limitation for the 
chordwise changes in angle of attack is because, if a chordwise change 
in angle of attack with a resulting load change occurs aft of any point 
of the wing trailing edge, the application of the two-dimensional 
vortex-sheet boundary conditions as used herein would give wing loading 
off the wing. Figures 3(a) and 3(b) indicate the limitations for two 
plan forms. The s panwise loading is independent of chordwise change of 
angle of attack in the shaded area of flgureB 3(a) and 3(b), and depends 
only on the final angle of attack attained at the line through the 
maximum upstream point of the trailing edge . The theory limitations 
require that the angle of attack of the wing surface, which is not shaded 
in figures 3(a) and 3(b), remain equal to the angle of attack attained 
at the line bordering the two surface regions. 

Figures 4(a) and 4(b) show swept wings approaching zero aspect ratio 
by modification of plan form in another way than that due to the Prandtl 
transformation. With the leading— edge sweep angle constant, the chord Is 
assumed to approach infinity. The two-dimensional vortex sheet can 
satisfy the boundary conditions for all such wings if the aspect ratio is 
made small in such a way that the leading— edge angle remains c ons tant «.nfl 
in the range -*/2 > < */ 1 2 - 

In summation, the plan forms for which the boundary conditions are 
satisfied by the two-dimensional vortex sheet for zero-aspect-ratio wirgs, 
or moderate-aspect-ratio wings at the speed of sound, include: 

1. All wings having an unswept trailing edge. (Spanwise loading 

is independent of chordwise angle-of-attack changes, depend- 
ing only on the angle of attack of the trailing edge, 
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provided the curvature at the trailing edge is small com- 
pared with the slope. The control— surface hinge line can 
have any position on the wing chord.) 

2. Those wings having swept trailing edge where every point of 
the trailing edge is at or "behind the upstream line of 
maximum wing span. (For the plan-form region at or ahead 
of the spanwise line through the maximum forward position 
of the trailing edge, the spanwise loading is independent 
of chordwise angle— of -at tack changes, depending only on the 
constant chordwise angle of attack of the region at or 
behind the spanwise line. The control— surface hinge line 
can have any position on the wing chord at or ahead of the 
spanwise line . ) 

The two-dimensional vortex theory employed herein is limited to 
unseparated flow, that is, flow that does not leave the wing surface 
along the chord. This requires that, for low aspect ratios, the aerody- 
namic angles be small. The extent that the aerodynamic results for 
control— surface hinge lines perpendicular to the air stream and blunt 
wings deviate from assumptions of the linearized— compress ible— flow 
equation at the speed of sound is not known. However, the aerodynamic 
results from the linearized theory for these cases are included in the 
present report since the results are independent of plan form for the 
plan— form range considered. The spanwise loading is also independent 
of the chordwise location of the control-surface hinge line for low-aspect— 
ratio control surfaces. 


Solution of Boundary Integral Equation 

Downwash due to the two-dimensional vortex sheet .— The boundary 
condition that, must be satisfied is that there be no flow through the 
surface of the airfoil. The downwash induced by the wing at, and normal 
to, the wing surface is thus equal to the upwash normal to the wing 
surface and experienced by the wing for a given angle— of-attack distribu- 
tion and free-stream velocity. 

The downwash due to a two-dimensional vortex sheet that extends to 
infinity fore and aft is given by 


v(y) 


_1 r b/,£ r«(y) 
^J-b/z 7-7 


ay 


( 1 ) 


With the assumption that the flow sticks to the wing plate, then the two- 
dimensional vortex sheet is at an angle of attack a and w(y) is the 
induced velocity normal to the wing plate. The upwash norm! to the wing 
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surface and experienced by the wing is equal to Y sin a(y) where 
a(y) is the spanwise angle-of-ettack distribution. The boundary condi- 
tion becomes 


w(y) = Y sin a(y) 
With the dimensionless relations 


( 2 ) 




y 

b/2 


equation (l) becomes with equation (2) and for small values of <x(y) 


a(n) 



(3) 


Fourier series solution of equation (3) Equation (3) gives the 
relation between the spanwise angle— of-attack distribution and the span- 
wise loading. Solution of the integral equation (3) for any spanwise 
distribution of angle of attack can give the vortex— strength distribu- 
tion and thus the spanwise loading. The solution of equation (3) can be 
obtained by a method similar to those of wing-section and lifting-line 
theories . 

The loading can be given by the trigonometric series 


G( cp) 


00 

^ a n sin ncp 
n=l 


w 


then 


0*(<P) 


00 

cos ncp 

n=l 


( 5 ) 


Equation (3) becomes with tj = cos cp and tj = cos cp v 
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a(<Pv> - r/ 


1 r * G*(<P)<1<P 


0 COB (p— COB 9v 




cos ncp 


cos cp-cos cpv 


dcp 


n=l 


The following Integral Is derived in reference 5: 


(6) 


P cos ncpdcp _ ii sin ncpy 
J cos cp-cos cp v sin cpy 


(7) 


Equations (6) and (7) give the solution of equation (3) as 


a(cpv) 


00 

-z 


na n sin nq> v 
sin cp v 


n=l 


( 8 ) 


Equation (8) relates the boundary conditions of the plan form with 
the unknown spanwlse load distribution. The boundary conditions are that 
the upwash experienced by the plan form across the span, given by 
a(cp v ), equals the downwash created by the circulation given by the series. 
The solution of equation (8) requires only that values of an be found 
that are independent of cpy . 

Equation (8) can be written as (omitting the subscript v since (pv 
is continuous) 


a(cp) sin cp 


00 



n B.JX sin ncp 


(9) 


which suggests the sine series of the Fourier series. The Fourier series 
for a function f(cp) is given by 


00 


f(cp) = 


n=l 


Bn sin ncp 
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where 



sin ncpdcp 


( 10 ) 


Let 


na^ = ^ 


( 11 ) 


then 



sin cp sin ntpdcp 


( 12 ) 


The Fourier series repeats the function f(cp) periodically for every 
2 rt interval. In the present problem we are interested only in the 
interval from 0 to it; then equation (12) becomes 


a n 



sin cp sin ncpdcp 


(13) 


With equation (13) , the values of the loading series given by 

equation (4) can be found for arbitrary angle— of-attack distribution. 


Evaluation of a n Coefficients 

The spanwise loading series coefficients can be evaluated for any 
reasonable function a(cp). Functions a(cp) that produce additional 
loading due to constant angle of attack, loading due to rolling, basic 
loading due to twist for zero total lift, and spanwise loading due to 
control-surface deflection will be considered. The a n coefficients 
that are determined can be inserted in equation (4) to obtain the 
spanwise loading distribution. 

Additional loading .- For additional loading the upwash experienced 
by the wing is constant spanwise and equal to Va, or a(cp) = a. Equa- 
tion ( 13 ) gives 


an. 



sin cp sin nplcp 
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a n = a for n=l 
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a n = 0 for n/l J ' 

Loading due to rolling .— For damping in roll the vertical velocity- 
experienced by the wing is proportional to the distance from the wing 
midsection and the rate of roll, or equal to -py = — V (pb/2V) tj, or 
a(cp) = — (pb/2V) cos cp. Equation (13) gives 


then 


"n 


-2(pb/2V) 

nit 



cos cp sin cp sin ncpdcp 


S = Z for n=e 1 

4 

> 

8^ = 0 for n/2 


(15) 


Basic loading .— Basic loading is due to twist and/or ‘ spanwise 
change in camber where, for zero aspect ratio, the camber is the slope 
of the midcamber line at the region of the wing trailing edge. As 
disti n guished from additional loading, the integral of the basic loading 
over the span is equal to zero, or Cl = 0. It will be shown that Cl 
is proportional to aj. hence for Cl = 0,- ai = 0. Equation (13) gives 
for n=l, and aa. = 0. 


£ 

it 



a(cp) 


sin 2 cp dtp 


0 


(16) 


Let 


a(cp) = arc + e(cp) 


(17) 


where a r is the angle of attack of the wing root chord for Cl = 0, 
and e(cp) is the wing twist relative to, the wing root chord. Equation 
(17) in (16) gives 


CGp — 



€(cp) 


sin 2 cp dep 


(18) 


With ocq- constant, equation (13) becomes with a^. = 0 for n>l 
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an 



sin cp sin nspdcp 


(19) 


The a n coefficients for basic loading can be determined from 
equation (19) for arbitrary s panvise tvist. The angle of attack of the 
wing root chord, for total lift equal to zero, is given by equation (l8). 

Basic loading due to linear tvist .— All the possible variations of 
tvist with span cannot be presented here. An example of the use of 
equations (l8) and ( 19 ) for the case of linear vashout or washin can 
readily be shown. For this case 

e(<p) = e 0 |cos cp ( 


where eo is a constant giving the total twist at the wing tips rela- 
tive to the root. A positive or negative value of e 0 gives washout 
or washin, respectively. 


Equation (18) gives 

_ _ ±lo r* /s 

a r “ « / 

J o 

Equation (19) gives 
for n >1 
-2e, 


cos cp sin 2 cp dtp = 


-4e c 

IF 


( 20 ) 


a n = 
then 



it/2 

cos <p sin cp sin ncpdcp — / cos cp sin cp sin ncp&cp 


■) 


a n = 0 for n=2 

— 4-e D sin (mt/2) 

VnA) for n >2 


(21) 


Two arbitrary control surfaces .— Two control surfaces that can 
(l) have any position on the wing span, (2) have unlike control-surface 
spans, and (3) be deflected differently can. be pictured as follows: 




62 

1 

It 1 

9 4 e 

3 «'/2 


PF 

Si 

v 

7 l i 

9 2 0 
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where 0 ’s are arbitrary wing-span positions of the control-surface 
edges; and Si and 62 are arbitrary control-surface deflections 
measured parallel to the plane of symmetry. 

Then in equation ( 13 ) 


a( cp) =0 from cp=0 to 0 X 
a(cp) =s Si from 0 i to 02 

a(cp) =0 from &z to 0 3 

a(cp) = 62 from 0 3 to 0 4 

a( cp) =0 from 0 4 to jf 


28i n 9 ? 28 2 r 6 * 

a n = EJT" / sin 9 sin + ST" / sin CP s3jl ^cpdcp 

J 6 X Jq s 


for n=l 




Sr 


02-01— ^ (sin 202-sin 20 x ) 




B 2 


0 4 — 0 3 — — (sin 20 4 — sin 203 ) 


} 


for n >1 


a„ = 


itn( n 2 — lj 


{ 6l 


cos 0i Bin n0i — cos 0 2 sin n0 2 


- > ( 22 ) 


n(sin 0i cos n0i-ein 02 cos n0 2 ) 


8 2 


cos 0 S sin nSg-cos 0 4 sin n0 4 — 


n(sin 0 S cos n0 3 — sin 0 4 cos n0 4 ) 


J 
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Symmetric spanvlse distribution of two control surfaces .— Two 
nn-ntrnl Rnrfftnas c\ -t k+.-H bntgifi RyTnrnfitvT f.s.lly on tEe wing flpan~nan be 
pictured as 


it 


I 


S 2 



02 



For this case 


©4 — IT — ■ 0i 
0 3 = It — 0 2 


then equation (22) gives 
for n-1 




a i 


5i+S 2 


02 r " ■ 0i — ^ (sin 20 2 - s in 20i) 


for n >1 


a n = 


-2( 81-62 cos nit) 
itn( n 2 — 1 ) 


cos 0i sin n0i— cos 0 2 sin. n0 2 — 


> 


(23) 


n(sin 0i cos n©i— sin 0 2 cos n0 2 ) J 

Inboard control surfaces .— Two inboard control surfaces of equal 
span can be pictured as 



0 
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for 

for 
an = 

span 


for 

for 


For this case 


02= ^ and equation (23) gives 


n=l 


a i 


( 5 1+82)/ < 

~~ 5 ~ ~~~ ( it— 20 3 .+S in 2 . 8 x ) 




n>l 

— 2(&x-82 cos mt) 
itn( n 2 — 1 ) 


\(2k) 


mt 

n cos — + cos 0 x sin n0x — n sin 0i cos n 0 


0 


Outboard control surfaces .— Two outboard control surfaces of equal 
can be pictured as 






31 


,1 

- 82 





It- 

it/2 

02 

0 



5i 


For this case 0i=O, and equation (23) gives 
n=l 


ax = 


Sx"^2 

"ST" 


(202-sin 20 2 ) 


n> 1 


(25) 


an = 


2(&x-52 cos nit) 
itnCn 2 — l) 


cos 02 sin n02~ii sin 02 cos n02 


)J 
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Aerodynamic Characteristics in Terms 
of the Coefficients 


In the foregoing sections the loading series coefficients were 
found for the several types of spanwise loading. In the following 
section aerodynamic characteristics will be found that are applicable 
to any of the loading series coefficient of the types of loadings con- 
sidered. 


The aerodynamic characteristics in terms of the a n coefficients 
are tabulated as follows: 

Loading distribution .— 

00 

> a n sin ncp (k) 

n=l 



Lift coefficient.— 


C L = A 



sin cp dtp = ^ a x 


(26) 


Rolling-moment coefficient .— 


C l 



sin cp cos cp dcp = 


( 27 ) 


Induced drag coefficient .— To find the induced drag, use can be 
made of the principles of Muni (reference 6), in which the total induced 
drag is found by considering the inclination of the lift vector at an 
angle equal to one-half the ratio of the downwash in the wake to the 
free— stream velocity. The downvaBh angle in the wake is given by 
equation (8). 


CDi = 



a(cp) sin cp dcp = 



n=l 


(28) 
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Semispanwise center of pressure for the right wing panel .— 



(29) 


Semispanwise center of pressure for the left wing panel.— 


n+1 


P G(cp) sin cp cos <p d<p na.z+6 V ^ ^ — 

Jft l—i n 2 


odd n 


^c.p.- 


r-> 


sin cp dcp 


00 

2sai+8 ^ 


even n 


(-1) 2 n a n 
n 2 — 1 


(30) 


With equations (4) and (26) through (30) the loading distribution, 
lift coefficient, rolling moment, induced drag, and spanwise center of 
pressure can be found for additional, rolling, basic, or control— surface 
loadings considered in the previous section. For convenience, the 
equations giving the aerodynamic characteristics due to the various type 
loadings are tabulated in table I. 


delations for the Various Types of Loadings 

Co mbina tion of the types of loadings .— The loading series coeffi- 
cients are additive for various types of loading as implied by equa- 
tion (13). The type of loading is determined by the angle— of-attack 
function, and summation of the angle— of— attack functions is equivalent 
to the su mma tions of the load ing series coefficients. Thus, all aero- 
dynamic characteristics that depend linearly on the a n coefficients 
are additive. As indicated by equations (4), (26), and (27), loading 
distribution, lift coefficient, anr? ro lling moment are additive; whereas 
induced drag and center of pressure (equations (28), (29), and (30)) are 
not. To find induced drag and center of pressure^ for summations of 
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several types of loadings, the a n coefficients can "be added and 
induced drag and center of pressure recomputed with equations (28), 

(29), and (30) with the total a^. values. 

Control-surface loading characteristics .— The slope of the loading 
distribution at any control— surface edge becomes infinite. From 
equation (4), 


G*(e) 


na a cos 0 


n=l 


With any of the control surfaces considered, G ’(0) contains a term 

of the type which series diverges and sums to infinity. Thus, 

mathematically the loading distribution is vertical at these points j 
however, in plots of loading distribution this vertical Blope is hardly 
discernible . 

From the discussion of the previous section, the loading distribu- 
tion for a given control-surface— angle distribution is additive with the 
loading for any other control— surface-angle distribution. As an example, 
this gives a simple relation between inboard and, outboard flaps (5 i= 6 2 ) 
for the same value of 0. 


G(cp) ,, , = 5 sin 9 — G(9 ) . , . 

' Y outboard ' inboard 

(°V) . * _ (°J£) 

\ A / outboard ^ \ A- / inboard 


(3D 


Closed function for loading series .- The spanwise loading distribu- 
tion due to control-surface deflection is somewhat unwieldy in the Beries 
expansion. The infinite series can be summed to a closed function 
involving only elementary functions. The summation is mathematically 
detailed in the appendix together with the closed functions of the 
distributions of control surfaces considered. 

Compressibility .— The equations of the present theory (tabulated in 
table I) can be multiplied through by 3 to give compressible parameters 
of aerodynamic coefficients and wing geometry. Then, subject to the 
limitations of the linearized-compressible— flow equation, the theory 
remains valid for low values of the parameter |3A throughout the speed 
range. Thus, at sonic speeds for which p is small, aspect ratio can 
have moderate values. 
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RESULTS 

Comparison With Other Theories 

The present theory is in accord with the results of reference 1 for 
additional loading due to wing angle of attack as evaluated in table I. 
The additional loading lift coefficient and induced drag are identical 
with the results of reference 3 for the triangular wing at the speed of 
sound for any finite aspect ratio. Loading due to rolling and the 
damping— irtr-r oil coefficient tabulated in table I as determined by the 
present theory are the same as that derived in reference 2 for the 
triangular wing of zero aspect ratio. 

Reference 1 indicates that, in the low— aspect-ratio range, the 
aerodynamic characteristics dependent on spanwise loading vary muc h less 
with aspect ratio than do aerodynamic characteristics dependent on 
chordwise loading. In the present report the values of n.nd pc^-p 

deviate, due to increasing values of PA, from the results of other 
theories and are given for several plan forms by comparing the present 
theory with the subsonic theory of references 7 R nd 8 and shown in 
figures 5 and 6, respectively. Figures 7 and 8 give the comparative 
ranges of PCt.^ and PCl as given by supersonic-theory values of 
references 9 and 10. 

The range of PA that the present theory approximates can be 
judged by the use of figures 5 through 8, and all span loading charac- 
teristics that are given by the present theory in addition to and 

Clp (e.g. , Clq and Cjg for control surfaces of low, compressible 
aspect ratio and full-^wing-chord control surfaces) are expected to be 
equally accurate in this PA range. For supersonic speeds the plan form 
influences more strongly the range of PA approximated by the present 
theory. As seen in figures 7 and 8 only the pointed wings in supersonic 
flow asymptotically approach the present theory, while blunt wings 
approach the present theory more abruptly. 


Figure Representation of Equations 

The span loading distribution due to inboard or outboard flap 
deflection with no differential in the flap angles is given by letting 
&i = S 2 = §f in equations ( 2k ) and (25) and use of equation (4). Since 
the poor convergence of the span loading series for control surfaces 
requires that n be taken to the order of 30 for the order of accuracy 
shown in the plots, the closed function of the loading distribution 
given in the appendix and the relation between inboard and outboard 
flaps given by equation (31) saves considerable computing time. The 
symmetric spanwise loading due to inboard and outboard flap deflection 
is given in figures 9(a) and 10(a), respectively, for several flap 
spans. As presented, 5 is measured parallel to the plane of symmetry. 



20 


MCA TN 2011 


For the limiting case of zero aspect ratio or at the speed of sound for 
moderate aspect ratio, G/Sf is constant for any fraction of wing-chord 
control surface, provided the control surface itself is effectively of 
low aspect ratio. 

The fact that spanwise loading is independent of the chordwise 
position of a control-surface hinge line is analogous to considering 
that the effective change of wing angle of attack to change of control- 
surface angle da/dS is unity for all fractions of the wing chord. 

Since da/dS is unity, the use of da/dS haB little significance for 
low-aspect-ratio control surfaces considered in the present report. 


The loading distribution due to inboard or outboard aileron deflec- 
tion with no differential between the ailerons is given by letting 
Si = -82 = S a in equations (24) and (25) and use of equation (4). The 
antisymmetric loading due to inboard., or outboard aileron deflection is 
given by figures 9(b) and 10(b), respectively, for several aileron spans. 
Presented also in figure 10 is the loading distribution for additional 
loading and that due to rolling . 

The induced drag coefficient of the wing due to inboard or outboard 
control surfaces is presented in figure 11 as a function of control- 
surface span. 

The lift coefficient due to symmetrically deflected inboard or 
outboard flaps is given by letting 5i = 62 = 5f in the respective 
equations of table I. Figure 12 gives the lift coefficient as a 
function of flap span and shows greater lift obtained with inboard 
flaps. 

The rolling-moment coefficient due to antisymmetrically deflected 
inboard or outboard ailerons is given by letting &i = — 82 = 8 a in the 
respective equations of table I. Figure 13 gives the rolling-moment 
coefficient as a function of aileron span and shows greater rolling 
moment obtained with outboard ailerons. 

Basic loading depends on the twist of the wing. Camber is equiva- 
lent to twist if the slope of the mean camber line at the region of the 
trailing edge is considered as the twist angle. For linear washout or 
washin, twist as given by e(cp ) = e 0 | cos cp|, the loading distribution 
given by the respective equation of table I is plotted in figure 14. The 
angle of twist at the wing tip given by- e 0 is positive for washout and 
negative for washin. 


CONCLUDING REMARKS 

The theory developed herein will satisfactorily predict the span 
loading of low— aspect-ratio wings that have all points of the trailing 
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edge at or behind the upstream line of maximum span. The same results 
apply to moderate-aspect— ratio wings at the speed of sound wherever the 
principles of the linear ized-compressible— flow theory are acceptable. 

Application of the theory shows that the spanwise loading, such as 
additional loading due to angle of attack, that due to rolling, basic 
loading due to twist, and that due to arbitrary control surfaces of low 
aspect ratio, is independent of plan form for the range of wings con- 
sidered. 

For the plan forms having unswept trailing edges, the spanwise 
loading is independent of the chordwise position of the control surface 
and depends only on the slope of the trailing edge. The total loading 
due to camber depends only on the slope of the trailing edge, provided 
the curvature at the trailing edge is s mall compared with the slope. 

For plan forms having swept trailing edge, the spanwise loading is 
independent of chordwise position of control surfaces, provided control- 
surface hinge line is at or ahead of a straight span line through the 
point of the maximum forward position of the trailing edge. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Aug. 15, 19^9. 


APPENDIX A 

CLOSED FUNCTION FOR CONTROL SUBFACE 
LOADING DISTRIBUTION 

The spanwise loading distribution of the loading due to control 
surfaces can be put into a closed function. The loading distribution 
is given as the series 


G(<p) = 


00 

^ a n sin ncp 
n=l 


where the an coefficients for the arbitrary control surface are given 
in equation (22). 
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Mathematics of the Series Summation 

For a specified point 0 at which a discontinuity in angle of 
attack occurs, the loading distribution is a function of two series 
given by 


G (<P) = (±) 


!{<«-! 


sin 20) sin cp + 


cos 


Z sin n0 
n(n= 

n=2 


sin ncp 
- 1 ) 


sin 0 


00 



cos n0 sin ncp 
n 2 ^-! 


} 


(Al) 


where the positive sign is taken for 0 measured at the left edge of 
the wing section of angle of attack 5, and-the negative sign for 0 
measured at the right edge of the wing section of angle of attack 5. 
In order to present G(<p) as a closed function, the indicated two 
summations of equation (Al) must be made. 


By the use of trigonometric identities the two infinite series of 
equation (Al) can be written, respectively, as 


and 


00 

1 V 1 cos n(6-<p) _ 1 

2 n(n2-l) 2 

n=2 


Z cos n(0-Kp) 
nCnili) 

n=2 


(A2) 

» 


00 00 

A. V flto MB+Sl _ 1 V sin n( 0-ff) 
2 n2-l 2 ns-l 

n=2 n=2 


(A3) 


Let 


(A4) 


z = e-^ = cos \|r + i sin if 
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Then expressions (A2) and (A3) are the real and imaginary parts of the 
respective summations given by 


00 


I 


z ± 


ntn 2 — 1) 


(A5) 


and 



n=2 


With the series. 


Zn(l- z) 


* -(z + 



(A 6) 


expression (A5) can be put into equivalent elementary functions by the 
following process: 


V 

n=2 



1 

z 



-Zn(l-z) 


1 dzdz 


3z 1 (1-s) 2 

T " 2 2z 


Zn (l-z) 


(A7) 


Similarly, expression (a 6) becomes 


00 



11=2 


i 

z 


/ 




n-1 


dz 


n=2 


—z Zn (l-z) 


dz 


1 , Z , 1— Z 2 / \ 

— + v + — ~ — Zn (1—z) 


2 i ’S + '2T 


(A8) 
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With equation (Ak) equations (A7) and (a 8) become, respectively. 


I ^L) ( 2 !n2 ! s±n 2 1 - l) eIn2 | 

n=2 


<* 


i( $ sin t + 2 sin 2 |[ cos 1 sin ^ 


(A9) 


and 


X n^j = ^ ( 2 + cos t) + sin \|r cos 1 sin ^ + 


n=2 


i (sin i|r) (j-j - 2 | sin | | ^) 


(A10) 


The s umma tions of expression (A2) are obtained by taking .the real, part 
of equation (A9) and are equal to 


— - (l-cos 9 cos cp)ln 


sin 


0-Kp 

~~T 


sin 


0-<P 

2 


1 

2 


sin 9 sin cp 2n 



+ 


3 

l 


sin 9 sin cp 


(All) 


The summations of expression (A3) are obtained by taking the imaginary 
part of equation (A10) and are equal to 


jj- cos 9 sin cp — sin 0 cos cp In 


£_ 

0-q> 

sin -=-• 


— cos 9 sin cp In 



(A12) 
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The quantity within the brackets of eqaution (Al) is equal to the 
product of cos 0 and expression (All) together with the product of 
—sin 0 and expression (A12) . The loading distribution becomes 


G(cp) 





0 sin cp + (cos cp — cos 0) Zn 

sln^ 



(A13) 


where, as before, the positive sign is tak^p. for Q measured at the left 
edge of the wing section of angle of attack S, and the negative sign for 
the right edge. Equation (A13) can be applied directly for the case of 
control surfaces considered in the text. 


Two arbitrary control surfaces .— 1 

sin^ 




G(cp) = £ 

(02-01) sin cp + (cos cp - cos 0 2 ) Zn 

2 

2 


(cos cp — cos 0 X ) Zn 


sin 


6 x 4<p 


sin 


6v-cp 


r 


2E. 

K 


(©4-0 3) sin cp + (cos cp — cos 0 4 ) in. 

03+9 


(cos cp — cos 0 3 ) in 


sin 


0^+2 


sin 


sin 


2 

63-cp 




(Aik) 


Symmetric spanwise distrib it ion of two control surfaces .— 

/c- ,c. \ r ' !_._ ®2+cp 

0(9) = - \ ^ (62-61) sin cp • ~ 


(cos cp — cos 0 2 ) Zn 


sin 


sin 


©2-CP 


(cos cp — cos 0 X ) Zn 


sin. ^ 


sin 


2 

— 


82 

it 


(cos cp + cos cp 2 ) Zn 


cos 


02+Cp | 


COS 


Qg-tp “(cos cp + cos cpi) Zn 


61+cp 


cos 2 


0i— cp 

cos — 

2 


- 


(A15) 
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Inboard control surfaces .— 


G(q>) = (^T^) (| “ 0 i) sln 9 


5i 

it 


cos 9 In 


1+ sin 9 
cos *9 


- (cos 9 — cos 0i) 2n 


sin 


sin 


01+9 

2 

na 


62 

It 


cos 9 In 


1-sin cp 
cos 9 


— (cos tp + cos 0i) Zn 


cos 


cos 


6,1+9 

2 

13 


(A16) 


Outboard control surfaces.— 


G( cp) = e 2 sin 9 + 


(cos 9— cos 0 2 ) In 


c 

-r (cos 9 + cos 0 2 ) In 



(A17) 
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TABLE I.- EQUATIONS FOR VARIOUS AERODYNAMIC CHARACTERISTICS 


BOIM- I BJKmvliS 

tion I loading distribution 





^ an sin cap 

For Cl ■ 0 


n«2 


(19) 


Or - — £ j «(q>) sin* <p<J* 

Jo 


«n-j ra ( «(<p) simp sin ncjxJcp 
Jo 

f 



Baaio loading due to linear twist 



m n~i 

k 'o V t— 1) 8 sin ap* 

Tor C L - 0 

(21) 

* ^ n(n=— V) 

n-3,5,7 


(where « 0 is positive for 

a. .-iSl 


washout, negative for waahln) 

** 3* 



(22) 

m 

^ an sin up 

5 (A14) 

n=l 


Loading due to two arbitrary control surfaces 8 


(ila 2 02 - i In SSi) j + ^ (Bin 0 - Bln® 5i) + 

(aln 25, - iln 26 a ) j ^ (Bln 5 S, - »ln“ fl a ) 


Loading due to syametrio distribution of two control surfaces 



§ (8i+6 a ) j^0ar0i - 



S (Si-»aK»ln s 9» - Bin 5 Sj) 


Loading due to inboard control surfaces 


(24) 

3 (A 1 S) 

M 

^ a n “In nip 
n«l 

^ (8i+6a)(»-fi9i + sin 29 x ) 

■g ( 81 — 6 a ) (1 — sin* 81 ) 

Loading due to outboard control surfaces 

(23) 

a (A17) 

m 

^ a n sin nip 
n«l 

4 (0i+S a )(2fl 3 - aln 25a) 
4 

~ (&i- 6 a) sin 8 0 a 


1 The span vise center of pressure for the right and left wing panels is given by equation [29} and ( 30 ) . 

2 &i and 0 + give the ving span pool tion (cos 8) of the outboard edges of the oontrol surfaces; 6 g and e s denote the inboard edges. 
®01osed functions for the epanvise loading distribution of oontrol surfaces are given in the appendix. 
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Figure I.— Vortex sheet of tow- aspect-ratio wing and corre- 
sponding approximate two-dimensional vortex sheet. 



Figure 2 Plan forms having un swept trailing edges approach- 
ing zero aspect ratio by Increasing longitudinal coordinates 
of plan forms . 
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Arbitrary control surfaces 



Figure 3 Plan forms having swept trailing edge showing shaded 
regions where chordwise variations of angle of attack are per- 
missible in the present analysis. 



Figure 4.- Plan forms approaching zero aspect ratio that main- 
tain a constant leading- edge sweep angle . 
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Figure 5.- Comparison of lift - curve - slope parameter, /BG La , 
with calculations of reference 7 . 



J3A 


Figure 6. - Comparison of damping -in - roll parameter, , 
with calculations of reference 8. 
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fiA 

Figure 7. -Comparison of lift -curve -slope para meter , ; fi Cl q , 
with calculations of reference 9. 



fiA 


Figure 8.— Comparison of damping- in-rott parameter, fiC lp , 
with calculations of reference 9. 





(b) Antisymmetric aileron deflection . 


Figure 9.- Wing semispan loading distribution, per radian, due to inboard con 
trot surfaces for various control -surface spans. 
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Figure //- Induced drag, per radian squared, due tQ symmetrically de- 
flected flaps or antisymmetrically deflected ailerons as a function 
of control- surface span. 



- Figure 12 -Comparative (iff, per- radian, due to symmetrically deflected 
inboard or outboard flaps as a function of flop span , 
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Figure 13.- Comparative rolling moment, per radian, due to antisymmetri- 
cally deflected inboard or outboard ailerons as a function of 
aileron span. 



Figure !4- Wing semispan loading distribution, per radian, for wings with 
symmetric linear twist, ± c„ = T ~ a r . 
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